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]_.  Introduction 

This  program  was  initiated  in  October  1976  when  ARO  commenced  funding 
of  our  studies  on  the  improper  ferroelectric  gadolinium  molybdate  Gd2(Mo04)3 
(or  GMO).  The  work  carried  out  on  GMO  under  the  present  program  can  be 
regarded  as  the  "winding  up"  of  a  major  study  previously  funded  by  AFOSR. 

Subsequently,  we  moved  into  studies  of  incommensurate  phase  transitions 
in  insulators,  which  produce  improper  ferroelectricity  when  the  associated 
sinusoidal  modulation  "locks  in"  by  adjusting  its  wavelength  to  become 
commensurate  with  the  underlying  lattice,  i.e.,  the  wavelength  becomes  a 
simple  multiple  of  the  basic  lattice  constant. 

We  began  by  studying  potass ium  selenate,  KgSeO^,  which,  at  that  time, 
was  the  only  system  known  to  exhibit  this  behavior.  As  other  isomorphous 
compounds  were  reported  to  exhibit  similar  behavior,  we  made  preliminary 
studies  of  these  systems,  but  the  bulk  of  our  effort  was  absorbed  by  ex¬ 
tensive  studies  on  KgSeO^  as  it  was  found  to  have  increasingly  rich  and 
complex  behavior.  Valuable  insight  into  this  was  provided  by  a  theoretical 
study  on  rubidium  calcium  trifluoride,  RbCaF-j,  which  provided  a  first- 
principles  understanding  of  the  structural  instability  in  this  simpler 
system.  From  this,  certain  general  principles  emerged  regarding  instabil¬ 
ities  in  systems  containing  molecular  groups  which  we  were  able  to  apply 
to  KgSeO^. 

Our  final  major  piece  of  work  during  this  funding  period  was  the 
development  of  a  new  model  for  the  behavior  of  GMO  which  appears  to  provide 
a  unified  explanation  for  a  wide  range  of  puzzling  anomalies,  some  reported 
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by  others  and  some  by  ourselves,  which  hitherto  have  defied  explanation.  In 
a  very  direct  sense,  this  is  the  true  conclusion  of  our  work  on  GMO,  since 
it  follows  very  directly  from  our  earlier  work  on  this  system  where  we  re¬ 
ported  the  first  observation  of  a  second  soft  mode,  apparently  associated 
with  the  transition  from  the  low- temperature  (improper  ferroelectric) 
phase  to  the  high- temperature  (paraelectric)  phase. 

2.  Nature  of  the  Studies 

Our  program  is  a  combined  experimental  and  theoretical  effort,  the 
former  employing  principally  inelastic  light-scattering  techniques,  and 
the  latter  principally  involving  detailed  lattice-dynamical  calculations 
for  complex  structures. 

In  order  to  provide  a  brief  but  adequate  overall  description  of  our 
work,  each  major  aspect  of  the  work  will  be  described  separately  by  an 
introductory  sunroary  followed  by  the  abstract(s)  of  our  major  papers  in 
that  area. 

2.1  Studies  on  GMO 

These  consist  of  our  earliest  ARO-sponsored  work  on  the  Raman  scat¬ 
tering  from  this  system,  in  which  we  identified  the  presence  of  two 
apparent  soft  modes  associated  with  the  phase  transition  from  the  ferro¬ 
electric  to  the  paraelectric  phase,  and  our  most  recent  work  in  collabo¬ 
ration  with  Or.  Jan  Petzelt  of  the  Czechoslovak  Acadenty  of  Sciences,  which 
advances  a  unified  theory  of  the  structural  instability  In  this  system. 

They  are  described  more  fully  in  the  following  two  abstracts: 
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«l  REVIEW  B  VOLUME  IS.  NUMBEB  7  Q|  3  * 1  OCTOBER  197S 

Aaomaloas  Auspbig  of  phonons  in  ferroelectric  Gd^MoO^ 

Q.  Kim*  and  P.  G.  Uliana 

MtkUn  Laboratory  of  Physics,  Unin rsity  of  Nebraska.  Lincoln.  Nebraska  (UN 
(EwM  21  April  1971) 

)lci  eauatr  of  the  A,  to—  tpectivru  of  ferroelectric  Gd^MoOJ,  ae  •  function  of  uniaxial  (treat  op  to 
,10*  |j  kbar  at  laap  if  tat  near  the  433JC  traaritioa  I-,  and  at  a  Amctioa  of  temperature  from  20  to 
)2)  K.  have  provided  aew  iaformatioa  (boat  the  anomalously  damped  tinea  peaking  at  about  30  and  75 
at  MO  K.  The  forocoaetaat  parameter  m,  of  the  75-cm'1  fine  a  independent  of  ttrem  whereat  its 

rn - r  dtewatee  with  iaereatrag  (treat  at  found  prcvioutly  for  the  50-cm'1  line.  No  other  linee 

yf  amoved  to  he  watt  dependent  in  tlrit  atrett  range.  The  (treat  dependence!  of  the  SO-  and  73-cm*1  linet 
yt  identical  An  aanlyrie  of  the  integrated  ktmtttirt  at  taro  ttrett  from  10  to  523  K  showed  ail  lines  except 
the  two  at  50  and  75  cm-1  to  approaimete  the  expected  Boon  Einetetn  temperature  dependence.  The  50-  and 
734m*1  hnet  damp  anomeloualy  on  heating;  their  integrated  intentitwe  go  through  maxima  near  320  K.  At 
15  K.  the  75-cm'1  hoc  tpfeti  into  a  doublet  with  II-  and  13-cm“‘  peaks,  timilar  to  the  known  splitting  of  the 
feem-1  fem  below  300  K  into  44.5-  and  51.5-cm*1  pcakt.  Only  two  hnet  below  T,  ate  allowed  by 
tjmmarf  to  cerwlete  with  foe  "toft"  doubly  degenerate  zone-boundary  phouont  of  the  perndectric  phase. 

Smte  it  them  eii  tiurreimti  none  of  the  ham  "eoftca**  with  incrcaring  temperature  or  ‘'harden’’ 

wah  incteetiag  ttrem,  the  idmritim  of  the  “toft"  phonone  in  foe  ferroelectric  phene  remain  undetermined. 


CENTRAL  PEAK-SOFT  MODE  COUPLING  IN  FERROELECTRIC  Gd^MoO^ 

J.  Petzelt  and  F.  Smutny 
Institute  of  Physics,  Czech.  Acad.  Scl. 

18200  Prague  8,  Na  Slovance  2,  Czechoslovakia 
V.  Katkanant,  F.  G.  Ullman  and  J.  R.  Hardy 
Department  of  Physics,  University  of  Nebraska 
Lincoln,  Nebraska  68588-0111 

and 

A.  A.  Volkov,  G.  V.  Kozlov  and  S.  P.  Lebedev 
Institute  of  General  Physics,  Acad.  Sci.  U.  S.  S.  R. 

117  942  Moscow,  Ul.  Vavilova  38,  U..S.  S.  R. 


iSVJ 


ABSTRACT 


Transmission  measurements  on  Gc^MoO^  In  the  5-50  cm”1  region  were 
performed  using  tuneable  backward-wave  oscillator  sources  (5-30  cm'1)  and 
a  Fourier  spectrometer  (30-50  cm'1).  The  resulting  dielectric  spectra  show 
an  additional  low-frequency  dispersion  which  was  fitted  to  a  standard  central 
peak  model.  Its  characteristic  relaxation  frequency  Is  ^20  cm'1  and  the 
coupling  between  the  soft  mode  and  central  mode  Increases  near  the  transition 
temperature.  This  model  also  accounts  very  well  for  the  weak  anomaly  In  the 
clamped  permittivity  ec  measured  at  63  MHz.  The  same  central  mode  was  used 
to  fit  earlier  Raman  soft-mode  spectra.  All  of  these  data  were  fitted  to  a 
three  coupled-mode  model  which  revealed  that  the  soft  mode  spectrum  consists 
of  two  strongly  coupled  bare  modes:  a  higher  frequency  mode  which  softens 
and  carries  the  entire  Raman  strength  and  a  lower  frequency  mode  which  Is 
hard  (59  cm'1)  and  Raman-Inactive.  Both  these  modes  are  also  coupled  to 
the  central  mode  and  this  coupling  Increases  sharply  near  the  transition. 

This  model  Is  consistent  wl  th  all  existing  experiments.  The  relatively 
large  width  of  the  central  mode  Indicates  its  intrinsic  nature  and  suggests 
partial  disorder  near  the  transition. 


2.2  Theoretical  Work  on  KgSeO^ 

The  Initial  object  of  this  work  was  to  develop  a  lattice-dynamical 
model  for  this  complex  structure  which  would  trace  the  origins  of  the 
Incommensurate  (modulation)  instability  to  a  balance  between  the  inter- 
ionic  forces.  This  was  successful  as  is  described  in  more  detail  in  the 
following  abstract: 


PHYSICAL  REVIEW  ■  VOLUME  21.  NUMBER  1  ((\AH-5)  I  JANUARY  1980 

Theoretical  lattice-dynamical  studies  of  the  incommensurate  phase 
transformation  in  KjSeO* 

M.  S.  Haque*  and  1.  R.  Hardy* 

Btt Mm  Labaratary  tf  Physics.  Uniwrsin  of  Nebraska.  Lincoln.  Nebraska  68SM 
(Received  ?  May  1979) 

W*  have  made  detailed  lattice-dynamical  undies  or  the  ori|in  of  the  incommensurate  phase 
transition  in  potassium  selenate.  These  were  hosed  on  a  ri*id-ion  model  with  Coulomb  interac¬ 
tions.  and  also  short-rinse  interactions  between  restricted  numbers  of  close  neighbors  The  as¬ 
sociated  short-ransc  force  constants  were  determined  using  the  static  equilibrium  conditions  for 
the  crystal  and  the  observed  Raman  frequencies.  The  calculations  were  made  for  the  room- 
temperature  structure  using  the  quasiharmonic  approximation  in  which  the  effects  of  varying 
temperature  were  simulated  by  varying  the  short-ransc  force  constants.  In  this  way  we  were 
abk  so  show  that  the  crystal  has  s  low-frequency  optic  branch  of  I,  symmetry  which  displays 
softening  and  instability  for  wave  vectors  Tin  the  vicinity  ofT-O.JT*  where  T*  is  the  first 
radptocal-kitice  vector  along  the  (1001  direction.  This  is  in  agreement  with  experimental  neu¬ 
tron  acattering  results  and  implies  a  transition  to  an  incommensurate  phase.  Full  sets  of  disper¬ 
sion  curves  at  130  K  for  the  (1001  direction  are  presented.  We  ako  present  comparisons  or  the 
— *  and  experimental  tow-frequency  Sj  and  !j  branches  at  130. 14S.  and  2S0  K.  We 
show,  by  decomposition  of  the  squares  of  the  normal-mode  frequencies  into  a  sum  of  Coulomb 
and  short-range  components,  that  the  balance  between  these  two  kinds  of  interatomic  forces  is 
very  delicate.  This  decomposition  also  shows  that  a  transition  to  an  incommensurate  phase  is 
very  likely  far  both  the  sekoatc  and  isomorphous  structures.  Our  model  also  enabks  us  to 
pndkt  that  unkxkl  stress  wiB  rapidly  depress  the  transition  kmperature;  a  result  also  in  accord 
with  experiment 


2.3.  Studies  of  Stress  and  Pressure  Effects 

As  the  last  abstract  explains,  we  found  that  the  balance  between  the 
Interlonic  forces  was  very  delicate,  suggesting  (as  preliminary  experimental 
work  had  demonstrated)  that  the  transition  temperature  could  be  markedly 


altered  by  the  application  of  both  uniaxial  stress  and  hydrostatic  pressure. 
This  led  to  a  series  of  experiments  designed  to  examine  the  stress  and 
pressure  dependence  of  both  the  incommensurate  and  commensurate  transition 
temperatures.  In  the  course  of  this  work,  two  collaborations  were  developed, 
one  with  Dr.  J.  D.  Axe,  at  Brookhaven  National  Laboratory,  in  which  neutron 
scattering  measurements  were  used  to  probe  the  effects  on  the  transition 
temperature  of  relatively  low  hydrostatic  pressures,  and  one  with  Dr.  6. 
Samara  at  Sandia  Laboratories,  in  which  dielectric  measurements  were  used, 
to  extend  the  range  of  pressures  to  much  higher  levels.  The  highest  pressure 
used  was  sufficient  to  suppress  the  commensurate  transition  entirely. 

This  work  is  described  in  the  following  abstracts,  along  with  our  own 
most  recent  studies  in  which  Raman  scattering  measurements  of  the  peak 
frequency  of  the  soft  mode  were  used  to  monitor  the  effect  of  uniaxial 
stress  on  the  incommensurate  transition  temperature. 


PHYSICAL  REVIEW  B  VOLUME  22.  NUMBER  I  JULY  19SO 

Effect  of  hydrostatic  pressure  on  the  incommensurate  phase  of  KgSeOa 

W.  Press 

Brookhaven  National  Laboratory.  Upton.  New  York  II97J 
and  Intinit  Jttr  Festkorperforschunp  det  K rrnfonchunttanlaitr. 

Jllhch.  West  Germany 

C.  F.  Mgjkrznk  and  J.  D.  Axe 
Brookbam  National  Laboratory.  Upton.  New  York  1197} 

J.  R.  Hardy,  N.  E.  Masse,  and  F.  O.  Ultman 

University  of  Nebratka-Lincoln.  Lincoln.  Nebraska  69598 
(Received  12  February  1910) 

A  neutron  scatterins  study  of  the  effect  of  hydrostatic  pressure  on  the  incommensurate  phase 
transformation  in  KjSeO,  was  performed  in  order  to  obtain  information  concerning  the  nature 
of  the  competing  forces  that  (ive  rise  to  a  structural  instability  at  a  general  wave  vector,  it  was 
found  that  an  increase  in  pressure  lowers  the  transition  temperature  ldT,/dF,  -  -6.5  K/kber) 
and  results  in  a  smaller  value  of  the  initial  incommensurate  wave  vector.  The  pressure  depen¬ 
dence  of  the  incommensurate  to  ferroelectric  phase  transition  was  also  investigated  (dTJdF, 
m—9  K/kbar).  The  latter  transition  displays  a  distinct  pressure  hysteresis. 
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VANISHING  OF  THE  PHASE  TRANSITIONS  AND  QUANTUM  EFFECTS  IN  K2Se04  AT  HIGH  PRESSURE 


G.  A.  SAMARA  . 

Sand  la  National  Laboratories!. Albuquerque,  New  Mexico  87185*  U.S.A. 
N.  E.  MASSA  and  F.  G.  ULLMAN 

University  of  Nebraska,  Lincoln,  Nebraska  68588,  U.S.A. 


Abatract-The  effects  of  pressure  on  the  dielectric  properties  and  phase  transition 
in  K2SeO«  were  investigated.  The  lnconmensurate-to-commensurate  ferroelectric 
transition  vanishes  at  0.73  GPa,  and  the  paraelectrlc-to- incommensurate  transition 
is  expected  to  vanish  at  ^  1.13  GPa.  These  large  pressure  effects  can  be  under¬ 
stood  in  terms  of  the  delicate  balance  between  Coulomb  and  short-range  interactions. 
The  results  show  the  increased  influence  of  quantum  effects  on  approaching  the 
dlsplacive  limit  (i.e.  Tc  -  0  K),  and  they  confirm  the  existence  of  a  y  -  2  regime 
in  the  temperature  dependence  of  the  susceptibility,  as  predicted  by  recent 
theories. 


Solid  State  Communication*,  Vo 1.42,  *0.3,  pp. 173-178.  1982. 
printed  in  Crest  Britain. 
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UNIAXIAL  S1USS  DEPENDENCE  OP  THE 
neotOKUBATE  PHASE  TRANSITION  IN  EjSeO^ 

N.  E.  Naaaa,  F.  C.  Oilman**,  and  J.  X.  Hardy 

Lam  Laboratory  of  Physics,  University  of  Nebraska,  Lincoln,  n«hra,fc»  68588-0111 

(locelved  16  November  1981  by  A.  A.  Maradudln) 

The  effect  of  uniaxial  stress  along  the  orthorhombic  "c"  axis  of  K*SaO, 

CrVltill  on  Pha  naalr  #vaa..aaa>  .  a  k.L. _ M _ a  ..  .  as  v2  4 


^  scattering  at  82  ,  93  ,  97  and  105  degK.  (The 
rl£  er*Mltlon  J »e  V-  «3K  *nd  the  lncomaeaaurate  transi¬ 

tion  is  la  the  range  Tj  -  123-130X,  at  aero  stress.)  The  change  in  Ti 
with  stress  is  estlaatsd  by  extrapolation  of  isobarlc  Curie-Heias 

•*  **  *9“r'  of  the  amplitude  mode  peak  frequency.  « J,  va.  T. 
a  “  'S'**  to  decrease  with  stress  along  the  "e"  axis  at  a  rata  of 
5  10  d*«X/kbar,  a  value  that  is  much  larger  than  the  -6.5 
tonaA  wlth  hydrostatic  pres  aura.  The  Curle-Walss  plots  show 
significant  sub-linear  deviations  at  thslr  low  temparature  ends.  Tbs 

t#  bBlow  300  hstsi  shove  300  bars,  the 

cryiuii  jritldid. 


2.4  Raman  Scattering  Studies 


Our  major  achievement  In  this  area  was  a  definitive  determination  of  all 
the  Raman  spectra  of  KgSeO^.  In  particular,  our  examination  of  the  high 
frequem  regl^  „  of  these  spectra  led  us  to  postulate  that  the  selenate 
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groups  are  orlentatlonally  disordered  above  the  incommensurate  transition 
(and  possibly  below).  Recent  Japanese  work  on  the  structure  of  isomorphous 
RbgZnCl^  has  demonstrated  such  disorder  to  exist  in  this  system  where  it  is 
quite  large,  and  thus  easier  to  detect. 

The  following  abstract  describes  our  KgSeO^  work: 
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Raman  scattering  from  KiSeO*  crystals  has  been  studied  in  the  (20— 800>-K  temperature 
mage.  Three  portions  of  the  spectrum  are  discussed:  defect-induced  scattering,  primarily 
below  100  cm-1,  the  external  mode  spectrum  below  200  cm-1,  and  the  internal  mode  spec¬ 
tra  in  two  regions,  300—500  and  800—950  cm-1.  The  temperature  dependence  of  the 
low-frequency,  defect-induced  scattering  has  been  correlated  (in  previous  studies)  with  the 
temperature  dependence  of  certain  nonzero-wave-vector  phonons  that  have  been  observed 
by  others  using  inelastic  neutron  scattering.  Close  to  the  incommensurate  transition  tem¬ 
perature,  129  K,  a  large  enhancement  and  line  narrowing  develops  below  10  ctn~';  no 
satisfactory  interpretation  of  this  effect  has  yet  emerged.  Above  129  K,  in  the  paraelectric 
phase,  the  multiplicity  of  the  external  modes  agrees  with  the  predictions  of  the  Raman- 
scattering  selection  rules.  The  internal  mode  spectra,  however,  consistently  contain  more 
lines  than  prescribed  by  the  selection  rules.  These  observed  internal  mode  spectra  would 
have  the  correct  multiplicity  if  the  crystal  did  not  possess  the  center  of  symmetry  of  the 
presumed  Pnam  space  group,  thus  suggesting  a  distortion  of  the  seienate  sublattice.  Small 
frequency  shifts  above  375  K  in  the  external  mode  spectra  are  similar  to  those  observed  in 
the  incommensurate  phase  and,  also,  differential  thermal  analyses  of  powders  and  crystals 
show  small,  diffuse,  reproducible  peaks  in  the  (373-473KK  range.  Doth  effects  may  be  as¬ 
sociated  with  glasslike  phase  changes  in  the  seienate  sublattice.  A  qualitative  model  of 
orientational  disorder  in  the  seienate  sublattice  is  offered  to  account  for  all  these  observe- 


2.5  Theoretical  work  on  RbCaF^ 


This  work,  whose  results  markedly  Influenced  our  understanding  of  the 
Raman  spectra  of  K2Se04,  was  a  collaborative  effort  with  Dr.  L.  L.  Boyer 
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of  the  Naval  Research  Laboratory  which  built  on  his  success  In  developing 
parameter-free  equations  of  state  for  the  alkali  halides.  It  represented  a 
major  extension  of  the  formalism  to  handle  the  more  complex  lattice  distor¬ 
tions  associated  with  the  cooperative  rotations  of  the  CaFg  octahedra  which 
are  responsible  for  the  structural  phase  transition  In  this  system.  RbCaF3 
was  chosen  as  the  trial  system  for  two  reasons:  It  was  the  system  about  which 
most  was  known  experimentally,  and,  more  Importantly,  good  free  Ion  wave 
functions  existed  for  all  the  constituents,  thus  making  possible  accurate 
Interionic  potential  calculations  for  all  pairs  of  Ions. 

The  following  abstract  describes  the  salient  conclusions: 
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We  have  made  a  flrst-principlea  study  of  the  structural  phaae  transition  at  Tt  - 193  K  in 
RbCaFj.  using  interionic  potentials  derived  by  the  Gordon-Kim  approach,  and  a  new  extension 
of  the  quasiharmonic  approximation  for  the  free  energy.  The  transition  is  caused  by  instability 
of  a  triply  degenerate  X-point  vibration  which  leads  to  a  coordinated  rotation  of  the  CaF,  oc¬ 
tahedra.  We  find  that,  as  the  lattice  contracts,  the  quasiharmonic  frequency  of  the  X-point  vi¬ 
brations  becomes  imaginary  at  approximately  1210  K:  Below  this  temperature  the  hoik  lattice 
energy,  as  a  function  of  CiF,  rotation,  has  a  double  minimum.  However,  the  quasiharmonic 
fret  energy  has  no  minimum  for  finite  rotations  until  T  £  125  K.  Thus  the  present  theory 
predicts  that  Tt  at  12J  K  (cf  Te  - 193  K,  experimental).  In  the  region  between  I2S  and  about 
1210  K  “neats"  of  modes  about  the  zone  edges  have  imaginary  quasiharmonic  frequencies.  By 
a  simple  extension  of  the  quasiharmonic  theory  their  contribution  to  the  free  energy  has  also 
been  included.  We  also  predict  that  the  melting  temperature  is  approximately  1350  K,  which 
agrees  very  well  with  the  measured  value  of  1312  K.  However,  the  predicted  thermal  expansion 
of  the  perovskite  phase  at  room  temperature  is  ~17%  tower  than  the  observed  value.  This  leads 
us  to  argue  that  the  good  agreement  between  theoretical  and  experimental  melting  temperatures 
is,  in  pert,  due  to  a  cancellation  between  neglected  anharmonic  effects  and  certain  deficiencies  in 
the  interionic  potentials.  We  also  find  that,  for  the  tetragonal  phase,  the  calculated  da  ratio  and 
rotation  angle  for  the  CaF*  octahedra  which  minimize  the  static  energy  are  in  good  agreement 
with  measured  values  at  low  temperature.  We  also  discuss  certain  more  general  implications  of 
the  present  work.  Specifically,  we  suggest  that  our  results  indies  te  that  it  may  be  more  natural 
to  regard  the  structural  phase  transition  u  arising  from  the  “unfreezing"  of  the  distortion  asso¬ 
ciated  with  the  lower-symmetry  phase.  Our  results  also  provide  a  natural  explanation  for  the  ap¬ 
parently  universal  tendency  of  transition  temperatures  for  zone-boundary  instabilities  to  be 
raised  by  hydrostatic  stress. 
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3.  Conclusions 

During  the  present  funding  period  we  have,  by  means  of  our  joint  ex¬ 
perimental  and  theoretical  studies,  proceeded  from  an  Initial  situation 
where  we  had  a  limited  understanding  of  the  structural  transition  In  one 
complex  system  (GMO),  based  on  limited  experimental  data  and  semi -empirical 
theory,  to  the  point  where  we  are  now  predicting  and  explaining  from  first 
principles  a  much  wider  range  of  observations  for  a  general  class  of  systems. 
We  now  have  the  ability  to  study  the  whole  range  of  structural  instabilities 
which  are  possible  when  displacement  patterns  of  non-zero  wavelengths  become 
unstable:  of  these,  the  transition  in  GMO  is  just  one  special  case;  in¬ 
commensurate  transitions,  as  typified  by  the  behavior  of  KgSeO^,  are  of  a 
much  more  general  nature  and  offer  the  possibility  of  a  rich  diversity  of  new 
physical  phenomena  whose  scientific  and  technological  Implications  have  yet 
to  be  explored. 


